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We propose a set of maximized string order parameters to describe the hidden topological order 
in the valence bond solid states of quantum integer spin-S chains. These optimized string order 
parameters involve spin- twist angles corresponding to Zs+\ rotations around z or x axes, suggesting 
a hidden Zs+i x Zs+i symmetry. Our results also suggest that a local triplet excitation in the 
valence bond solid states carries a Zs+i topological charge measured by these maximized string 
order parameters. 
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Determining order parameters is one of the most im- 
portant issues in the study of strongly correlated sys- 
tems. It is the basis of Landau theory in describing con- 
tinuous phase transitions between different phases with 
spontaneous symmetry breaking. However, some novel 
phenomena or phases, for example, the fractional quan- 
tum Hall states and the Haldane gap phenomena in quan- 
tum integer-spin chains, are not amenable to such a de- 
scription. In particular, to find an order parameter in a 
system with strong frustrations or quantum fluctuations 
is challenging and highly nontrivial jj, |2j ■ 

For quantum integer-spin antiferromagnetic Heisen- 
berg chains, Haldane predicted that there is a finite exci- 
tation gap above the ground state [3[ . To understand this 
intriguing conjecture, Affleck, Kennedy, Lieb, and Tasaki 
(AKLT) proposed a set of exactly solvable quantum spin 
models whose ground states are the valence bond solid 
(VBS) states |4j]. In the 5" = 1 VBS state, there is a hid- 
den antiferromagnetic order which can be characterized 
by string order parameters (SOPs) 5] . Corresponding to 
this nonlocal string order, Kennedy and Tasaki showed 
that there is a hidden Z2 x Z2 symmetry [6|. To reveal 
this symmetry clearly, they introduced a nonlocal unitary 
transformation to convert the VBS state to a "diluted" 
ferromagnetic state and showed that the ground state is 
fourfold degenerate if this symnretry is broken in an open 
chain. 

For quantum integer-spin chains with S* > 1, it is ex- 
tremely difficult to find the corresponding nonlocal uni- 
tary transformation to convert the VBS state into a long- 
range ordered state. So far no one has succeeded along 
this direction. However, this is not the only approach 
for revealing a hidden topological symmetry. This hid- 
den symmetry can be also identified from the non-local 
string order parameters of the ground state. So the ex- 
tended SOPs with a twist angle 9 = n/S were proposed 
to account for the hidden order in the VBS states [71, Isj. 
Actually, this kind of extended SOPs does not give a 
comprehensive description to the Haldane phase in these 



systems [9|. Recently, we have extended the 5 = 1 VBS 
(which is SO{3) invariant) to the SO{2n + 1) symmetric 
system [10|. We find that these S0{2n + 1) symmet- 
ric matrix product states are the exact ground states of 
Hamiltonians of spin S = n with either S0{2n -|- 1) or 
SU{2) symmetry. Furthermore, we show that this system 
possesses a hidden topological (Z2 x Z2)" symmetry and 
the ground state can be characterized by n-set of hidden 
antiferromagnetic orders. This study suggests that the 
topological order is generally more complicated than the 
0-twisted SOP as previously studied in the VBS states. 
It motivates us to consider how to find proper SOPs to 
characterize the hidden topological order in the higher 
integer-spin VBS states [ll|. 

In this paper, we will propose a systematical method 
to find the optimized SOPs in the VBS states. By gen- 
eralizing the S' = 1 den Nijs-Rommelse SOP to arbi- 
trary integer spin-5 VBS states, the corresponding hid- 
den long-range order can be most comprehensively man- 
ifested by the maximal string correlation functions be- 
tween two polynomial spin operators 



^" = E «" (s?) 



(a = x,y,z), 



(1) 



n=0 



where a„ are the coefficients which maximize the string 
correlation function. We find that these maximized SOPs 
suggest a hidden Zs+i x Zs+i symmetry, in consistent 
with the degeneracy of the ground states in an open chain 
case. Moreover, our optimized SOPs measure a Zs+i 
topological charge carried by a localized triplet excitation 
in the VBS state. 

The AKLT Hamiltonian is defined by 0, ^ 



2S 



-ffAKLT^E E Kjrjii,l + 1), 



(2) 



i j=s+i 



where Kj > and Pj(j, z -I- 1) is to project two spins at 
i and i + I onto the subspace of the total spin J. This 



Hamiltonian can be written as a polynomial of SU{2) in- 
variant nearest-neighbor spin exchange interactions. The 
VBS state is the ground state of this Hamiltonian. 

In the Schwinger boson representation, the spin op- 
erators are expressed by S^ = ajbi, S~ = bjai, Sf = 
{ajai — h\bi)/2 with a local constraint a\ai + b]bi = 25. 
The VBS state in a length-i periodic chain is then given 

by 



|VBS) = n(al6l+i 



t„t 



b a 



i+ll 



(3) 



i=l 



In an open chain, the VBS states have effectively two 
spin-5/2 edge states at the two ends of the chain. In 
the thermodynamic limit, L —^ oo, the VBS states with 
different edge states are asymptotically orthogonal to 
each other, giving rise to (5*+ l)^-fold degenerate ground 
states. 

The VBS state can be also represented in a matrix 
product form [8l.ll3l.ll4|. By using the binomial theorem. 



it can be shown that 

s 
|VBS) = ^ Tr(B["il---B[™^l)|mi---mL), (4) 

{mi} = -S 

where i?[™l is a (S* + 1) x (S* + 1) matrix defined by its 
matrix elements as 

5M(p,5) = {~lf-P+^^{S + my.{S-m)\ 




m,q—pi \^) 



where 1 < p, q < S + 1. 

The correlation function of the VBS state in the matrix 
product form can be evaluated using the transfer matrix 
method [8|, [ij]. To do this, let us first introduce the 

following {S + 1)^ X {S + 1)^ transfer matrix: 



where P is a local operator acting on a single site, and 
B denotes the complex conjugate of B. For the iden- 
tity operator P = I, Gp is a Hermitian transfer matrix 
G* = Em (-^'™' <8)Pt"l). With these definitions, it is 
straightforward to show that the two-point spin correla- 
tion function can be expressed as 

<,b,b^)-iim^ Tr(G)i ' ^> 



where 



Gs = ^ TO (^["1 ® B[" 



In the long distance limit, the two-point spin correlation 
function always decays exponentially with the distance 
between the sites i and j, 



{sts^) 



exp 



\j-^\ 
e 



(8) 



where ^ = l/ln(l + 2/5) is the correlation length [12| . 

In order to describe the hidden topological order in 
the VBS states, let us introduce the following generalized 
string correlation function [8|: 



i-i 



OliO)^ lim ((^n^n^'''''"^")' 



(9) 



where y4" is defined by Eq. ([1]). The value of oq in A" 
can be fixed by demanding the expectation value of A" 

to be zero, i.e., ao = - En=i «« {{SfT)- 

Since the VBS state is spin SU(2) rotational invariant, 
we need only to evaluate the z-component of the string 
correlation function 0\{9). Based on the transfer matrix 
technique, it can be shown that 

Tr[(G)^-J+'-iG^t(Gop-''-iGA] 



0^^(61) = lim lim 

|j— i|— >oc L— »oc 



Tr(G)^ 



(10) 

where Go, G^t, and Ga are obtained by replacing op- 
erator P in Eq. dSl) by exp(i6'5^), (A^)^ exp(i6'5^), and 
A^ , respectively. We emphasize that only G and Go are 
Hermitian, while G^t and Ga are not. 

In the limit of |j — i| ^ oo and L ^ oo, O\{0) is deter- 
mined purely by the largest eigenvalues and eigenvectors 
of G and Gq- Wc find that Eq. PH)) can be simplified as 



oue) = 



1 



X2 



(•^maxl Ga |Aniax)| 



(11) 



where |Amax) and jA^ax) are the eigenvectors correspond- 
ing to the same largest eigenvalue of G and Go , respec- 
tively. 

When 6 = 0, 0\{6) becomes the ordinary two-point 
correlation function. It vanishes because lA^j^^) = jAmax) 
for 6' = and the vector Ga |Amax) is orthogonal to 
|Amax)- If we fix the form of the operator A^' and tune 
the spin- twist angle 6, 0\{9) will reach its maximum 
only when the vector Ga |Amax) is parallel to the vector 
l-^max)- This means that in order to maximize the SOP, 
lA^a^) is orthogonal to |Amax), i-e., (A^axl^max) = 0. For 
the VBS state, the eigenvectors |Amax) and lAj^^x) can be 
readily calculated, we find that the orthogonality condi- 
tion can be satisfied if and only if the following equation 
is satisfied 



l + e" 



128 



e'^' = 0. 



(12) 



Thus, the spin-twist angle corresponding to the maximal 
SOP is determined by 

2n7r 



5 + r 



(n=l,---,5). 



(13) 



Importantly, these spin-twist angles correspond to a dis- 
crete symmetry group Zs+i- Because of the spin SU(2) 
rotational symmetry, a similar maximal SOP can be 
found in terms of the x-component operators O^ {6) . The 
spin-twist angles corresponding to the maximal 0^{9) 
are also given by Eq. p^ . 

In the conventional Landau theory, it is well-known 
that the maximized order parameters fully characterize 
the symmetry of the low-temperature ordered phases, 
and these order parameters can be also used to describe 
the possible phase transitions as decreasing the tempera- 
ture from the high temperature disordered phases. When 
we generalize the similar rule to the present case, the 
maximal SOPs fully describe the hidden long-range or- 
der of the spin- 5 VBS states, suggesting that there ex- 
ists a hidden Zs+i x Zs+i symmetry. In a finite open 
chain, this hidden Zs+i x Zs+i symmetry is broken in 
the ground states, leading to the spin-|^ edge states with 
{S + l)^-fold degeneracy. We would like to emphasize 
that these spin-|^ edge states are dictated by the under- 
lying low-energy effective field theory, i.e., 0(3) nonlinear 
sigma model plus a topological term jl9|. Moreover, the 
hidden Zs+i x Zs+i symmetry also incorporates the pre- 
vious results given by Oshikawa [7|. For an odd integer 
S chain, the Zs+i group has always a Z2 subgroup with 
9 = 7:. However, for an even integer S chain, a twist 
angle ^ = tt is always absent. This result explains why 
the hidden Z2 x Z2 symmetry described by the den Nijs- 
Rommclsc SOP is broken in the odd-S* VBS states but 
not in the even-S* ones [7]. 

Let us consider the hidden topological symmetries in 
the first few cases of the spin-S VBS states. For the S = 1 
VBS state, we have A^ = S^. The maximal eigenvalues 
of the transfer matrices G and Go are equal to Amax = 3, 
and their eigenvectors are given by 



|A 



max/ = 4^(1 1)^, 



V2 
1 



|ASax) = ^{e^' ly , 

respectively. With a simple calculation, it can be shown 
that 0%{9) = (4/9)sin^(6'/2), from which its maximal 
value corresponds to 9 = tt. This result is fully con- 
sistent with the orthogonality condition. According to 
the nonlocal unitary transformation, the corresponding 
VBS state does exhibits a ^2 x Z2 topological symmetry. 
Therefore, it is a reliable method of using the maximal 
SOPs to reveal the hidden symmetry. 

Next let us consider the 5 = 2 VBS state, the maximal 
eigenvalues of the transfer matrices G and Go is Amax = 
40 and the corresponding eigenvectors are given by 

lAmax) - -^(1 1 1)^, 
|A2ax) = ^{e''' e'« l)^ 



When e = 27r/3 or 6* = 47r/3, the SOP dH) is maximal. 
From the condition that Ga lAmax) is parallel to |Amax)j 
we can find the corresponding optimal operators display- 
ing the maximal string correlation function as the form 



^] - ^3 



± 



7 



m' - 2] , 



(14) 



for 9 = 2tt/3 and 9 = 47r/3, respectively. The maximal 
value of the SOPs are O^ (9) = 3. In this nontrivial 
case, it is the linear combination of spin operator Sj and 
its spin quadrupole operator (Sj)"^ — 2 that exhibits the 
maximal hidden string order instead of the spin operator 
itself. The hidden topological symmetry in the S = 2 
VBS state corresponds to the Z3 x Z3 discrete symmetry. 
For the S = 3 case, the spin- twist angles for the max- 
imal SOPs are determined hy 9 = 7r/2, tt, 37r/2 from the 
orthogonality condition. Similar to the previous proce- 
dure, it is straightforward to show that the corresponding 
combinations of the spin operators are given by 



A' 



S^ 



± ^ mr ■ 



for 9 



67 

7r/2 or 37r/2, and 

89V2 



4] 



7 
67 



{S. 



2\3 



A', 



67 



S' 



14 



(sn 



z^3 



(15) 



(16) 



for 9 = TT. The maximal value of the SOPs is 0\ (9) = 
2592/4489, and the corresponding hidden symmetry is 

Z4 X Z4. 

Accordingly, the above discussions can be readily ex- 
tended to any higher integer spin-S' VBS states. In ad- 
dition to the VBS states, it can be also generalized to 
reveal the possible hidden topological order of arbitrary 
matrix product states. Thus, this is a systematic ap- 
proach to analyze the hidden topological symmetry of 
matrix product states. 

In fact, the topological property of the ground state 
can be also understood from its elementary excitations. 
For the AKLT model, to create an elementary excita- 
tion is to insert a triplet defect in the VBS ground state. 
This kind of excitation is called crackion [l6|. Under 
the Kennedy- Tasaki unitary transformation for 5* = 1, 
a crackion is a kink linking two ferromagnetic ordered 
ground states from left to right. In higher-S' systems, al- 
though the Kennedy- Tasaki unitary transformation does 
not transfer the VBS state to a ferromagnetic one, the 
kink or soliton-like feature of crackion can still be re- 
vealed by studying the string correlation function in the 
presence of a triplet defect [8| . An example of a crackion 
in the S = 2 VBS state is depicted in Fig. [1] 

With the Schwinger boson representation, the wave 
function of a localized crackion between k and fc + 1 
sites can be constructed by replacing a singlet (a^6[,^_j — 
with one of the triplet operators T° 



bWk+i 



Ti 



,t„t 



a'k^'k+i^ ^° = 4^1+1 + &I4+1' T 1 = blbl^^. 



The corresponding matrix product wave functions for 
these excitations are 

\^l)= Y^ Tr(B['"il---Ci™^-l---B['"^l)|mi---TOL), 

(17) 
where Ca is defined by 



l„ 



C^^\p.q) =^{-lf~''^/{S + m)\{S-m)\ 
C'-Iip.^)] 

For one crackion in the string, it can be shown that 
the SOP in the thermodynamic hmit is given by 



|j-«l 






(A„ 



"Vmaxl (j'^ |Amax) 

(18) 
where 0\{9) is the maximized string order parame- 
ter without crackions. When replacing operator P by 
exp(j0S'^) and / in the foUowing crackion transfer ma- 
trix: 

G?, = ^ (m'l P \m) (Cl"'l ® C^l) , (19) 



one obtains Gq 



definitions, 

■max I "-^O 



and G", respectively. From these 
one can show that the ratio between 
(ASaxl^o lA^a^) and (A^axl G" |Amax) is equal to 
exp(za6') (a = ±1,0) for the three kinds of crackions. 
If there are a few diluted crackions between i and j, the 
above result can be extended to 



j-i 



hm ((Ant-Q^,:.S^^. 



k—i k—i 

(20) 
where the phase factor exp{i9 J2k=i ^k) counts the total 
number of crackions between i and j. 

For the 5 = 1 VBS state, we have 6* = tt and the SOP 
alternates its sign according to the parity of J2k=i'^k, 
which can be interpreted as a Z2 topological charge [8|. 
If we consider the SOP in the x-direction, a crackion 
with a = also flips the sign. In a general spin-S* case, 
the spin-twist angle of the maximized SOP is given by 
e = 2Trn/iS+l) (n = 1, . . . , 5). Our optimized SOP sug- 
gests that the crackion carries a Zs+i topological charge, 
in consistent with the hidden Zs+i x Zs+i symmetry ar- 
gument. 

In summary, we have shown that the hidden order 
in the VBS states can be characterized by the general- 
ized den Nijs-Rommelse-type SOPs. The maximization 
of these SOPs automatically leads to spin-twist angles 



FIG. 1: (Color online) The schematic of a crackion excitation 
in a spin-2 VBS state. The solid lines represent the valence 
bond singlets. A dot represents a spin-1/2 Schwinger boson. 
The four bosons enclosed by each dashed circle form a spin-2 
state. The arrow represents a crackion of local triplet defect. 



corresponding to Zs+i rotations around 2; or a; axes, sug- 
gesting the existence of a hidden Zs+i x Zs+i symmetry. 
In the presence of the crackion excitation, the maximized 
SOPs are shown to exhibit a Zs+i topological charge. 
Recently, it was shown that the den Nijs-Rommelse SOP 
is an effective measure of the localizable entanglement in 
the S = 1 VBS state [17|. We believe that our maxi- 
mized SOPs provide a natural extension of the den Nijs- 
Rommelse SOP and can be used to explore multipartite 
entanglement properties of arbitrary VBS states. 

We thank D. H. Lee and X. G. Wen for stimulating 
discussions. Support from NSF-China and the National 
Program for Basic Research of MOST-China is acknowl- 
edged. 



[i; 
[2: 
[3; 
[4; 

[5; 

[6; 

[7: 
[s; 

[9; 
[lo; 

[11 

[12 
[13 

[14 



* Electronic address: lgmzhang@tsinghua. edu.ciil 
S. Furukawa, G. Misguich, and M. Oshikawa, Phys. Rev. 
Lett. 96, 047211 (2006). 

D. Perez-Garcia, M. M. Wolf, M. Sanz, F. Verstraete, 
and J. I. Cirac, Phys. Rev. Lett. 100, 167202 (2008). 
F. D. M. Haldane, Phys. Lett. 93A, 464 (1983); Phys. 
Rev. Lett. 50, 1153 (1983). 

I. Affleck, T. Kennedy, E. H. Lieb and H. Tasaki, Phys. 
Rev. Lett. 59, 799 (1987); Commun. Math. Phys. 115, 
477 (1988). 

M. den Nijs and K. Rommelse, Phys. Rev. B 40, 4709 
(1989). 

T. Kennedy and H. Tasaki, Phys. Rev. B 45, 304 (1992); 
Commun. Math. Phys. 147, 431 (1992). 
M. Oshikawa, J. Phys.: Condens. Matter 4, 7469 (1992). 
K. Totsuka and M. Suzuki, J. Phys.: Condens. Matter 7, 
1639 (1995). 

U. SchoUwock, O. Golinelli, and T. Jolicoeur, Phys. Rev. 
B 54, 4038 (1996). 

H. H. Tu, G. M. Zhang, and T. Xiang, J. Phys. A 41, 
415201 (2008); Phys. Rev. B 78, 094404 (2008). 
S. Michalakis and B. Nachtergaele, Phys. Rev. Lett. 97, 
140601 (2006); C. Hadley, Phys. Rev. Lett. 100, 177202 
(2008). 

D. P. Arovas, A. Auerbach, and F. D. M. Haldane, Phys. 
Rev. Lett. 60, 531 (1988). 

M. Fannes, B. Nachtergaele, and R. F. Werner, Europhys. 
Lett. 10, 633 (1989); Commun. Math. Phys. 144, 443 
(1992). 

A. Kliimper, A. Schadschneider, and J. Zittartz, J. Phys. 
A 24, L955 (1991); Z. Phys. B: Condens. Matter 87, 281 
(1992). 



[15] T. K. Ng, Phys. Rev. B 50, 555 (1994). Phys. Rev. Lett. 92, 087201 (2004); L. Campos Venuti 

[16] G. Path and J. Solyom, J. Phys.: Condens. Matter 5, and M. Roncaglia, Phys. Rev. Lett. 94, 207207 (2005). 

8983 (1993). 
[17] F. Verstraete, M. A. Martin-Delgado, and J. I. Cirac, 



